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The computation of higher order processes very often involves a large number of 
diagrams. In addition, it is in general not possible to solve the occurring integrals 
explicitly and expansions in small quantities have to be performed. This makes it 
necessary to automate the calculations as much as possible. A program package will 
be described which generates automatically the Feynman diagrams, manipulates 
the expressions in the desired way and performs the computation. As a physical 
application 0(aa s ) corrections to the decay rate of the Z boson into bottom quarks 
are discussed. 



1 Introduction and motivation 

The impressive experimental precision reached so far mainly at LEP, SLC 
and TEVATRON has made it necessary to increase also the effort from the 
theoretical side. A crucial role plays thereby the computation of multi-loop 
diagrams in order to evaluate quantum corrections to the different observables. 

Higher order corrections are mostly accompanied with a large number of 
diagrams. Very often it is a tedious job to do the bookkeeping and not to forget 
some relevant contributions. Especially when the number exceeds the order 
of a few hundred the generation should be passed to the computer. A further 
major task is the computation of the integrals. In most cases an exact solution 
is by far not possible and one has to rely on approximations. One possibility 
it to compute the integrals (at least partly) numerically. Another promising 
attempt is the use of asymptotic expansions which is applicable as soon as a 
certain hierarchy exist between the mass scales involved in the process. Well- 
defined prescriptions provide rules which specify the actions on the individual 
diagrams. In general each diagram generates several subgraphs which further 
increase the complexity of the calculation. Thus it is desired to automate the 
asymptotic expansion procedures. Of course, also the very computation of 
the single terms needs to be done by the computer as the size of intermediate 
expression become rather large. 
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In the next section a possible solution to the problems addressed above is 
presented by means of the package GEFICOM. In Section [3] its application to the 
computation of 0(aa s ) corrections to the decay of the Z boson into bottom 
quarks is discussed. 



2 GEFICOM: A package for generation and computation of Feynman 
diagrams 

The automation of the computation of Feynman diagrams can be divided into 
different steps. In this section we will discuss them on the example of GEFICOMEl. 
A flowchart is pictured in Fig. |l|. 

The user has to provide a few simple files specifying, e.g., the process, the 
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Figure 1: The structure of GEFICOM. 
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particle content and the allowed vertices. Then, in a first step, the graphs con- 
tributing to the considered process have to be generated where the amplitudes 
at least need to contain information about the involved particles and the mo- 
mentaJlowing through each propagator. Inside GEFICDM the FORTRAN program 
QGRAFd is used which has the advantage of being quite fast. Moreover different 
output formats are available one of which contains all the desired information 
on the diagrams. 

The Feynman diagrams are classified according to their topology which 
becomes especially important during the calculation of the integrals. This in- 
formation is, however, not provided by QGRAF and has to be determined from 
its output. For a human being this would be very easy. For a computer, on 
the other side, this is a quite non-trivial task and the momentum distribution 
has to be used in a clever way in order to extract the topology. This is done 
by FindTop which is part of a Mathematica program essentially inserting the 
Feynman rules and transforming the QGRAF output into FORMS notation. Ad- 
ditionally administrative files are generated which rule the computation, sum 
up the diagrams etc. One finally ends up with a huge database containing all 
relevant files. 

The very computation of the diagrams has to be. initiated by the user. It 
is based on two FORM packages, MINCERQ and MATADH, being able to deal with 
massless two-point, respectively, massive bubble diagrams up to three loops. 
Before, however, the amplitudes are passed to these packages there is the 
possibility to apply asymptotic expansions namely the hard-mass and large- 
momentum procedure providing consistent prescriptions on how to get expan- 
sions for large internal masses, respectively, large external momenta. They 
work on a diagram-by-diagram basis and make well-defined rules available 
which determine the subgraphs to be extracted from the original diagrams. 
The subgraphs have to be expanded in their small quantities before any mo- 
mentum integration is performed. Thus significant simplifications are obtained 
and very often the integrals can be written as products of lower order ones at 
the price of — sometimes significantly — increasing their number. The au- 
tomation of these procedures has been performed in the packages LMPQ and 
EXPlI. They apply the asymptotic-expansion procedures and express the initial 
diagram as products of single-scale integrals which can be treated by MINCER 
and MAT AD. 

A very powerful tool for the computation of the diagrams is implemented 
into MINCER and MATAD, the so-called integration-by-parts techniaufi. Up to 
now it has been systematically applied up to the three- loop level BEM. The 
general idea is quite simple: The integration-by-parts method exploits the fact 
that the momentum integrals in D dimensions are finite. Thus the surface 
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terms are zero: 



/d^/(P,...) = o. (i) 

The function / in general is a product of scalar propagators belonging to the 
considered loop. Performing the derivatives in Eq. ([!]) explicitly leads to a sum 
of several terms which has to vanish identically. Combining different equations 
of this type, which are obtained by varying f(p, . . .) and choosing different 
loop momenta, leads to relations — so-called recurrence relations — where 
a complicated integral is expressed as a sum of simpler ones. The repeated 
application of recurrence relations finally leads to a sum simple integrals, which 
can be reduced to integrals of lower loop order, and only a few, so-called, master 
integrals which require a hard calculation. 

At the three-loop level Jjhe described strategy has first been applied to 
massless twopppint integrals lI and has later on been extended to massive bub- 
ble diagramsE2l. Also the corresponding master integrals are availableE3. These 
two sets of recurrence relations were implemented into MINCER and MATAD, re- 
spectively, and we will see in the next section that they provide when combined 
with the asymptotic expansions mentioned above a very powerful tool with a 
large field of application. 

3 Application 

The properties of the Z boson have been measured to a very high accuracy 
- sometimes even at the permille level. Also the hadronic width of the Z 
boson, Thad, is known to an accuracy of roughly 0.1%. An important contri- 
bution to Thad is provided by the partial rate into bottom quarks as already 
at one-loop order the top quark enters as a virtual particle giving rise to en- 
hanced corrections proportiopaLxp GfM ( 2 . The full electroweak corrections axe 
known since quite some iimell^'E-a and also the leading terms of 0(a s GFM^)i3 
and 0{aa s lnM t 2 /M^)li3 are available since already a few years. Actually it 
turned out that the logarithmic contribution is relatively large as compared to 
the term quadratic in M t although \n(M? / M^) is by far less enhanced than 
M t 2 /M^. Thus the question arises whether the M t -independent coefficient 
and the power-suppressed terms in the expansion lead to a significant change 
in the numerical prediction. 

In the approach chosen intj the Z boson propagator was considered whose 
imaginary part directly leads to the decay width. The diagrams contributing at 
the two-loop level are shown in Fig. 0. A gluon has to be added in all possible 
ways in order to obtain the graphs leading to corrections of order aa s . 
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Figure 2: Diagrams contributing to Y(Z — > bb) and containing a virtual top quark. Thin lines 
correspond to bottom quarks, thick lines to top quarks, dotted lines to Goldstonc bosons 
and inner wavy lines represent W bosons. 

For vanishing bottom quark mass up to four scales may appear inside a 
single diagram: M 2 , M^, ^w^w an< ^ Q 2 where q is the external momentum to 
be identified with M§ at the end. £w is the electroweak gauge parameter. In 
order to be able to apply the hard-mass procedure a certain hierarchy has to be 
fixed. The assumption M 2 ^> (M^, q 2 — is clearly justified. Furthermore 
the diagram of Fig. ||(/) makes it necessary to choose ^> q 2 in order to 
avoid contributions to the imaginary part arising from cutting lines involving 
W bosons. At first sight this looks seemingly inadequate. However, a closer 
look to the corresponding diagrams shows that actually (2Mw) 2 , respectively, 
(M t + M w ) 2 has to be compared with q 2 which justifies the above inequality. 
This is true, since the real emission of a W boson is always accompanied by 
another W boson or a top quark. The notation ;§> q 2 is only formal, 
telling the hard-mass procedure which subgraphs have to be selected. 

Concerning the fourth scale, ^wM^, there is some freedom for the choice 
at which place in the inequality chain it is inserted. Only q 2 3> ^w^w is 
not allowed as then again unwanted imaginary parts would be generated. The 
other three possibilities are allowed and must lead to identical final results as 
the dependence drops out at the very end. In intermediate steps, however, 
the expressions which have to be evaluated are different. Thus the different 
choices provide quite strong checks on the calculation. 

At two-loop level seven diagrams have to be considered. Their calculation 
would still be feasible by hand. At three loops, however, 69 diagrams contribute 
and a computation by hand is very painful especially if higher order terms in 
the 1/Mt expansion are considered. The hard-mass procedure applied to the 
69 initial three-loop diagrams results in 234 subdiagrams which have to be 
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expanded in their small quantities. InliU the package EXP written in Fortran 90 
was used which automatically applies the hard-mass procedure. It can be called 
form GEFICOM setting appropriate options in the initial files to be provided by 
the user. 

Since we are interested in the virtual effect of the top quark we consider in 
the following the difference of the partial widths of the Z boson into bottom 
and down quarks: 

sr w d = 5r o,w _ sr o,w > (2) 



where #rj is the contribution from the diagrams involving a W boson to 
the partial decay rate T{Z — > 5T^ W is obtained from the diagrams in 

Fig. U and the corresponding ones at three-loop order. The zero indicates that 
next to the vertex diagrams no additional counterterms had been introduced 
as they would drop out in the difference. Note that besides the pole parts also 
the dependence on £w drops out in Eq. (Q). The final result reads in numerical 
formO: 
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r°lq - 0.50 -0.07+^ 



1.16 + 0.13 



(3) 



with r° = N c Mza/(12sgCg), sg — sin9w, &w being the weak mixing angle, 
c 2 , = 1 — and Ltw = ln(Mf /M^). Mt is the on-shell top mass. The numbers 
after the second equality sign correspond to successively increasing orders in 
1/M 2 , where the brackets collect the corresponding constant, lnM t and, if 
present, In 2 M t terms. The numbers after the third equality sign represent 
the leading Mf term and the sum of the subleading ones where the 0(a) and 
0(aa s ) results are displayed separately. 

One observes that for the realistic values M z = 91.91 GeV, M t = 175 GeV 
and Sg — 0.223, which were used in Eq. (||), the constant at next-to-leading 
order dominates over the In M 2 term known before. It is remarkable that at 
one-loop level the corrections arising from the 1/M 2 terms are of similar size 
than the one from next-to-leading order, however, the signs are different. The 
higher order corrections in 1/M t are smaller, which means that effectively only 
the leading M 2 term remains. 

Proceeding to two loops the situation is similar: Starting at C(l/M 2 ) the 
sign is opposite as compared to the leading terms and a large cancellation 
takes place. Here, the 1/M t 4 term is still comparable with the 1/M? contribu- 
tion. The 1/Mf term, however, is considerably smaller which suggests that the 
presented terms should provide a reasonable approximation to the full result. 
Comparison of this expansion of the one-loop terms to the exact result ofli-3 
shows agreement up to 0.01% which also gives quite some confidence in the 
0{aa s ) contribution. 

For a = 1/129 and a s (M z ) = 0.120 we get: 
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(-5.69 - 0.79 + 0.50 + 0.06) McV = -5.92 MeV. 



(4) 



The first two numbers in Eq. (||) correspond to the 0(a), the second two to 
the 0(aa s ) corrections. Each of these contributions is again separated into 
the Mf terms and the sum of the subleading ones. 
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